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ABSTRACT 



We quantify an evolutionary channel for single sdB stars based on mergers 
of binaries containing a red giant star and a lower mass main sequence or brown 
dwarf companion in our Galaxy. Population synthesis calculations that follow 
mergers during the common envelope phase of evolution of such systems reveal a 
population of rapidly rotating horizontal branch stars with a distribution of core 
masses between 0.32 M — 0.7 M Q that is strongly peaked between 0.47 M — 
0.54 M Q . The high rotation rates in these stars are a natural consequence of 
the orbital angular momentum deposition during the merger and the subsequent 
stellar contraction of the merged object from the tip of the red giant branch. We 
suggest that centrifugally enhanced mass loss facilitated by the rapid rotation 
of these stars may lead to the formation of single sdB stars for some of these 
objects. 

Subject headings: stars: horizontal-branch — stars: rotation — stars: subdwarfs 



1. Introduction 



Subdwarf B (sdB) stars are core He-burning st ars that haye an extremely thin (< 



0.02 M rn) hydrogen-rich envelope (for recent reviews, see lHebeiil2008l ; ICatelanll2007;lMoni Biden et al. 



2007 ) . Observations indicate that 40 - 70% of sdB stars in the field are in binaries (IMaxted et al 



2001 
2006). 
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Han et al.l (120021 . 120031 ) identified and performed extensive population synthesis calcu- 
lations of two channels for the formation of sdB stars in binaries in the field: (1) one or two 
phases of unstable Roche lobe (RL) overflow leading to common envelope (CE) evolution 
and (2) one or two phases of stable RL overflow, and one channel for the formation of sin- 
gle sdB stars in the field: the merger of two He white dwarfs (WDs). Their model of the 
population of sdB stars in binaries has been very successful in explaining several features of 
the observed field population, such as the prevalence of systems with P < 10 days, and this 
model is generally well accepted. Their model of the single sdB star population in the field 
pr edicts a relativel y flat distribution of masses ranging from ~ 0.4 M — 0.7 M (see Fig. 12 
in lHan et al.ll2003l ). This distribution differs from the mass distrib ution expected from sin- 



gle star horizontal branch (HB) models (e.g., iDorman et al.lll993l ). which predict a narrow 



distribution of sdB star masses strongly peaked at ~ 0.47 M Q . Observationa lly, masses for 



seven single field sdB stars have been determined using asteroseismology (see iRandall et al. 
20071 and references therein). Six are very close to 0.47 M Q and one is 0.39 M Q . However, it 



is not evident whether these masses are representative of the single sdB star population as 
a whole, since only ~ 5% of sdB stars pulsate. Further, reliable mass determinations using 
asteroseismology depend on correctly identifying the pulsational modes of the star, which 
can be quite difficult. 

The formation of single sdB stars also has been modeled using channels that do not 
require the merger of two WDs, but rather involve a single star with a phase of enhanced 
mass loss on the red giant branch (RGB). Various suggestions for the cause of this enhanced 
mass loss include He mixing driven by internal rotati on (Sweigar t 1997), spin up of the 
primary via interactions with a planet in a close orbit (jSokerl Il998l ) . and the existence of a 
sub-population of ZAMS stars with a high He-abundance (e.g., D'Antona et al. 2005). In all 
of these proposed causes, the amount of mass loss on the RGB must be fine-tuned in order 
to yield the very small envelope masses observed in sdB stars. However, it is not evident 
how such fine-tuning is incorporated naturally into the physical mechanisms that have been 
suggested to enhance the mass loss. 

In this Letter, we quantify a channel for the formation of single sdB stars in the field of 
our Galaxy: the evolution of objects resulting from the merger of stellar or sub-stellar objects 
with RGB stars during a phase of CE evolution initiated by unstable RL overflow or by a 
tidal instability (hereafter, referred to as the "CE merger channel"). We carry out detailed 
population synthesis calculations of these mergers and predict the resulting population of HB 
stars at the present epoch, some of which may have very small envelope masses and could be 
observed as single sdB stars. Although our treatment was developed independently, the idea 
that single sdB stars might result from the merger o f an RGB star and a low mass companion 
during CE evolution was suggested by lSokerl (119981 ) and explored further by lSoker fc Harpaz 
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f l2000l . l2007f ). 



Method 



We use the Monte Carlo population synthesis code developed by Politano (IPolitano 



19961 ; IPolitano fc Weilerl 120071 ). For the present study, two major mo difications have been 
incorporated: 1) the analytic fits that had been used previously (see iPolitanol Il996l ) have 
been replaced by numerical tables of 116 up-to-date stellar models ranging in mass from 0.5 
to 10.0 M in increments of 0.1 M and from 10.5 to 20 M Q in increments of 0.5 M Q and 2) 
tidal effects which act to synchronize the rotation of the primary with the orbit. 



The updated stellar models were ca 



culated using a version of the binary stellar evolution 



code, STARS, developed by Egglet on (|Eggletonlll97ll . Il972l ; lEggleton &: Kiseleva-Eggleton 



20021 ) and updated as described in iPols et al.l (119951 ). The treatment of convective mixing, 
convective overshooting, the helium flash and our definit i on of the core mass and of the 
envelope binding energyj are given in Ivan der Sluys et al.l ( 20061 ). Convective overshooting 



on the main sequence is taken into account for stars with M > 1.2 M Q and helium is ignited 
degenerately in stars with M < 2.0 M Q . The initial composition of our model stars is similar 
to solar composition: X = 0.70, Y = 0.28 and Z = 0.02. Mass loss via stellar winds 



is incorporated in the models using a R eimers-type pre s cripti on (iReimersi Il975l ) with an 
effective 77 of 0.2 and the prescription by Ide Jager et al.l (119881 ). which is negligible except 
for the most massive stars in our grid. 

The second modification to the code is that synchronism of the primary's rotation with 
the orbit is now assumed when the primary is on the RGB or the asymptotic giant branch 
(AGB). For these cases, after each time step, the new total angular momentum (AM) of the 
system is calculated, accounting for losses due to stellar winds. Synchronism is then enforced 
by redistributing this total AM such that the rotational angular velocity of the primary is 
equal to the orbital angular velocity of the binary. In all cases, we assume a circular orbit 
and rigid body rotation for the primary and we neglect the rotational AM of the secondary. 
We also assume that the rotational AM of the primary on the ZAMS is negligibly small. 

In our population synthesis calculations, we begin with 10 7 zero-age main sequence 
(ZAMS) binari es. We assume t hat th ese binaries are distributed over primary mass, M Pi n, 



according to a iMiller fc Scald (119791 ) IMF, over orbital period uniformly in log Pq (lAbt 



lr rhe gravitational potential and internal thermal energies of the envelope are included in the binding 
energy, but not the ionization/recombination energy. 
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19831 ) and over mass ratios uniformly in q (i.e., q(qn)dqo = ldqn), where go = M R n/M p n 
and M s n is the mass of the secondary star (jDuquennoy fc Mayorl Il99ll ; iMazeh et al.l 11992 



Goldberg et al.ll2003l ). Once the primary mass and mass ratio have been generated for a 
given binary, the secondary mass is determined by their product: M s0 = qoM Pi o. We adopt 
a minimum primary mass of 0.95 M Q , a maximum primary mass of 10 M and a minimum 
secondary mass of 0.013 M in our calculations. For secondaries with masses less than 0.5 M Q , 
including substella r secondaries, we use detailed stellar models from the Lyon group (e.g., 
Baraffe et al.l 120031 and refs. therein). For secondaries with masses greater than or equal 
to 0.5 M Q , we adopt the same stellar models used for the primaries. We assume that the 
duration of the CE phase is neg ligible compared to the other time scales (e.g. 



Iben fc Livio 



19931 ; iTaam fc Sandquistl 120001 ). that the age of the Galaxy is 10 yrs, and that the star 



formation rate throughout the Galaxy's history has remained constant. 

To model the population of mergers between giant primaries and less massive compan- 
ions, we employ the following scenario. A given ZAMS binary is evolved until the primary 
reaches the base of the RGB, at which point we assume synchronism between the primary's 
rotation and the orbit. As the primary ascends the RGB, CE evolution can be initiated 
either via unstable RL overflow or a tidal inst ability. We adopt t he criterion for unstable 
mass transfer via RL overflow given by eq. 57 in lHurley et al.l (120021 ) . To determine if merger 
occurs within the C E, simple energy considerati ons are used to relate the pre- and post-CE 
orbital separations (ITutukov &: Yunglesonlll979l ) according to a standard prescription (e.g., 
Willems et al.ll2005l ). The binding energy of the primary's envelope at the onset of the CE 
phase is determined directly from the stellar models. The efficiency at which orbital energy 
is transferred to the CE, parameterized by oce, embodies a major uncertainty in this sim- 
plified prescription. We have chosen «ce = 1 as our standard model and we investigate the 
dependence of our results on a CE . Merger is assumed to have occurred if the radius of the 
secondary is larger than its RL radius at the end of the CE phase. In addition, we consider 
mergers that result from a tidal instability. If the moment of inertia of the primary exceeds 
one-third of the moment of inertia of the binary before the primary fills its RL on the RGB 
or AGB, there will no longer be sufficient AM within the orbit to keep the envelope rotat- 
ing synchronously (jDarwmlll879l ) and the orbit will continue to shrink until the envelope is 
ejected or a merger occurs. 

Available detailed hydrodynamical models of stellar mergers have focused on high- 
velocity collisions between stars and are more appropriate for mergers that occur as a result 
of stellar interactions within dense stellar environments, such as at the center of clusters 



[e.g., iSills et al.lll997l . l200ll ; lLombardi et al.ll2002l ). In the absence of detailed models that 
are applicable to stellar mergers during CE evolution, we use a very simplified model for the 
merger process. The key assumptions in our model are listed below. 
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1. In most cases, the merger of the entire secondary with the primary leads to an object 
spinning at several times its break-up rotational velocity, given by Vhr = (GM/Rfl 2 , where 
M and R are the object's mass and radius, respectively. Consequently, we assume that a 
rapid mass loss phase occurs when the rotational velocity of the object exceeds a critical 
rotational velocity, v CT i t = |i>br, which corresponds to centrifugal forces contributing about 
10% of the force against gravity. 

2. We assume that the primary assimilates only as much mass from the secondary as 
is required such that it rotates at this critical velocity and we assume that the remaining 
secondary mass is expelled from the system. 

3. Although the merged object is not likely to be in a state of thermal equilibrium 
initially, we assume for defmiteness that the radius of the merged object corresponds to that 
of a star in the model grid for the same core mass and total mass after it achieves thermal 
equilibrium. 

4. In cases where choosing a model with the same core mass and total mass would result 
in an unphysical evolutionary regression (such as an RGB star reverting to a Hertzsprung gap 
star as a result of the merger), we assume that the radius of the merged object corresponds 
to a star in the model grid with the same total mass, but with the smallest core mass within 
the same evolutionary state as prior to merger. 

The subsequent evolution of the merged object, including its rotational velocity, v TOt , 
is then followed until the present epoch. If t> rot > t> cr it at any time, we assume that AM 
is removed via mass loss until the object is rotating at exactly f cr j t . We make the same 
assumptions regarding the choice of stellar model following mass loss as we did following 
mass accretion during merger. If the merged object completes its evolution as an AGB star 
before the present epoch is reached, we assume the merged object will be a white dwarf at 
the present epoch. 

3. Results 

In Figure 1, we show the two-dimensional distribution of the rotational velocities (as a 
fraction of v crit ) and the envelope masses (M env ) for our predicted population of stars on the 
HB at the present epoch. The distribution of rotational velocities rises slowly from 0.08t> cr it 
to ~ 0.35f cr it, remains fairly uniform up to 0.95f cr i t , and then rises rapidly, with 28% of 
the population rotating at v CI ft. The distribution of envelope masses peaks at ~0.35M Q and 
decreases with increasing envelope mass, following the IMF for the primaries. Approximately 
24% of the population has M env < 0.50 M and there are a small number of objects with 
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very low envelope masses, ranging from 0.27 M down to 0.07 M Q . 

In Figure 2, we show the distribution of core masses for our predicted population of HB 
stars at the present epoch (solid curve). We find that the core masses range from 0.32 M 
to ~ 0.7 M and the distribution is strongly peaked between 0.47 and 0.54 M Q , with 78% 
of the population having core masses in this range. For comparison, the dashed curve in 
Figure 2 shows the predicted distribution of core masses in a present-day population of HB 
stars that results from normal single star evolution. The two distributions are rather similar, 
although the dashed curve is slightly less peaked than the solid curve, with 67% of the dashed 
population having core masses between 0.47 and 0.54 M Q . Also, the number of present-day 
HB stars from our merger channel is about an order of magnitude smaller than the number 
predicted from normal single star evolution and the HB stars from our merger channel are 
rotating 10 - 30 times faster on the average than those from normal single star evolution. 

All of the HB stars in Fig. 1 are the result of a merger that occurred when the primary 
was on the RGB. During the merger, we have assumed that the primary ceases to assimilate 
mass from the secondary once the rotational velocity of the merged object reaches a critical 
rotational velocity, v cv - lt = | t^ r - Thus, by assumption, the merged object begins its evolution 
to the present epoch rotating at t^ crit . As the object evolves up the RGB, the increase in 
radius causes the moment of inertia of the envelope to increase and the rotational velocity 
to decrease correspondingly. The rotational velocity continues to decrease until the object 
reaches the tip of the RGB and He is ignited in the core. If He is ignited degenerately, the 
object's radius decreases by a factor of 5 - 15 and the moment of inertia of the envelope 
decreases dramatically. If He is ignited non-degenerately, the object's radius and moment 
of inertia still decrease, but much less dramatically. In the majority of cases, the decrease 
in the moment of inertia causes the object's rotational velocity to exceed v CI ft. Should this 
occur, we artificially remove mass from the envelope (at the specific AM of envelope's radius) 
before taking the next evolutionary time step, until the envelope is just rotating at f C rit- 

In 

approximately three-quarters of the population, the growth of the object's radius during 
core He burning is sufficient to eventually reduce the rotational velocity of the envelope to 
sub-critical by the present epoch. However, in the remaining one-quarter of the population, 
the envelope continues to rotate at v cr i t at the present epoch and rotationally enhanced mass 
loss may continue. 

We have investigated the sensitivity of these results to our assumed choices for a CE and 
g{%) by selecting two other choices for each: « CE = 0.5 and 0.1, and g(qo) = Qo an d Qo°' 9 - 
Selected results for these cases and for our standard model are given in Table 1. We find 
that the various features of the population described above do not depend significantly on 
the assumed choice for either acE or g (g ). The distribution of core masses remains strongly 
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peaked, with the fraction between 0.47 M Q and 0.54 M varying only slightly from 78% to 
84%. A significant fraction of the population still rotates at w crit , with this fraction varying 
from 29% to 39%, and has envelope masses less than 0.5 M , with this fraction varying from 
24% to 29%. In all models, a small number of systems have very low envelope masses, down 
to 0.07 - 0.08 M . The choices for «ce and g(qo) have a stronger effect on the total number 
of stars in the population, N, as indicated by the last column in Table 1. Compared with 
our standard model, N increases by as much as 50% for the «ce = 0.1 model and decreases 
by 56% for the g(qo) = q^ ' 9 model. 

We have also investigated the effect of increasing the assumed value of f cr i t on our 
population of HB stars. As shown in Table 1, choosing v cr i t = «br still results in a significant 
fraction of the population rotating critically and in a sharply peaked distribution of core 
masses. However, the fraction of the population with M cnv < 0.5 M is drastically reduced, 
resulting in no systems with M env < 0.35 M . As this choice for v cr i t is rather extreme, our 
results suggest that a critical velocity corresponding to a fraction of the break-up velocity is 
more appropriate. 



4. Discussion 

In our model population of HB stars, a small fraction (< 0.1%) have envelopes with 
masses less than 0.1 M (see Fig. 1). Approximately one-half of these systems have rotational 
velocities equal to v CT i t , suggesting that these stars will continue to experience enhanced mass 
loss. Follow-up evolutionary calculations of core-He burning stars with an initial envelope 
mass of 0.1 M (using the same STARS stellar evolution code with an artificially large wind 
to simulate rapid mass loss) suggest that at some point (~ 0.07 M ) the star will begin to 
contract upon further mass loss and move blueward along the HB. It is conceivable, therefore, 
that these critically-rotating, very low envelope mass HB stars in our population may lose 
sufficient envelope material to become sdB stars. 

These potential sdB systems constitute a very small fraction of our population. Upon 
examination of Fig. 1, this fraction would increase substantially if the envelope masses in 
our population were systematically reduced by ~ 0.2 M . Such a reduction might be accom- 
plished by assuming a higher rate of mass loss due to stellar winds on the RGB. We have 
adopted a standard Reimers-like prescription for mass loss in our stellar models, but it has 
been suggested that a rapidly rotating envelope may lead to enhanced mass loss on the RGB 
(e.g., ISweigartl 119971 ). In addition, we find that our assumption of synchronization on the 
RGB results in small number of primaries being spun up to t; cr i t prior to contacting their RL. 
This number is increased significantly if we choose ^ crit < | v hr . For example, if v crit = 0.1 tv, 
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two-thirds of the primaries are spun up to v CV i t prior to contacting their RL. However, we 
are not able to properly model either possibility in our population synthesis calculations, 
since a self-consistent treatment of mass loss in such cases would require stellar evolutionary 
models in which mass loss via winds is coupled to the rotation and this is beyond the current 
capabilities of our code. 

Reliable quantitative predictions concerning the absolute number of single sdB stars 
produced from our CE merger channel are not provided here, since our treatment of the 
merger process is too crude and is severely limited by the lack of models of stellar mergers 
relevant to CE evolution. Rather, our intent in this Letter has been to investigate a possible 
evolutionary channel for the formation of single sdB stars and to motivate its further study. 
In comparison with previously suggested single sdB star channels, this channel has certain 
advantages: (1) the AM necessary to spin up the envelope prior to the HB is e xplicitly and 
naturally provided by the orbital AM of the companion (see also ISokerl Il998t ) ; (2) the CE 
phase may be initiated at any point along the RGB, not just near the tip, since mass loss 
from the envelope also occurs on the HB as a result of super-critical rotation caused by the 
star's contraction following He ignition in the core; and (3) the predicted spectrum of core 
masses is consistent with that expected in HB stars. 

Further, our population synthesis calculations are, to the best of our knowledge, the first 
to explore in any substantive manner a possible observational counterpart to mergers result- 
ing from CE ev olution. It is generally accepted that merger s are likely to occur during CE 



evolution (e.g.. Ilben fc Livid Il993l ; iTaam fc Sandquistl 120001) and with a frequency roughly 



comparable to that of the surviving post-CE binaries (e.g.. iPolitano fc Weilerl 120071 ). Thus, 
there may exist a significant, but as yet unrecognized, population of objects in our Galaxy 
that are the result of a merger during CE evolution. Assuming this is true, identification of 
possible observational counterparts of these systems is likely to stimulate the development 
of models necessary to study them. 
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Table 1. Selected results for various input parameters 



Model 


Fraction of population with 


iV/iVstand 


0.47 M < M c 
< 0.54 M 


V TOt — 

v crit 


M env < 
0.5 M 


standard 


0.78 


0.28 


0.24 


1.00 


a C E = 0.5 


0.81 


0.32 


0.25 


1.19 


«CE = 0.1 


0.84 


0.39 


0.29 


1.50 


9(Qq) = 1o 


0.79 


0.29 


0.24 


1.01 




0.78 


0.31 


0.24 


0.44 




0.66 


0.21 


0.06 


1.10 
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Fig. 1. — Unnormalized two-dimensional distribution of envelope masses and rotational 
velocities in our model population of present-day HB stars from the CE merger channel. 
The rotational velocities are given as a fraction of the assumed critical rotational velocity for 
mass loss, v crit (see text). The total number of HB stars in our present-day model population, 
N, is 13,100. 
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Fig. 2. — Unnormalized distribution of core masses in our model population of present-day 
HB stars that results from the CE merger channel (solid curve). The distribution of core 
masses in a model population of present-day HB stars that results from standard single star 
evolution is also shown for comparison (dashed curve). The dashed population was calculated 
using the same stellar evolution models, IMF, star formation rate and total number of ZAMS 
stars as in our mergers population. The total number of HB stars shown in the dashed curve 
has been normalized to the same number as in the solid curve to facilitate comparison; in 
actuality iV S i ng = 8.5 iV m erg- 



